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Abstract 

The eigenvalues of an arbitrary quaternionic matrix have a joint probability 
distribution function first derived by Ginibre. We derive the j.p.d. for the 
weakly non-Hermitan version of this problem and then show that there exists 
a mapping of this system onto a fermionic field theory. This mapping is used 
to integrate over the positions of the eigenvalues and obtain eigenvalue density 
as well as all higher correlation functions for both the strongly and weakly 
non-Hermitian cases. 

I. INTRODUCTION 

Several ensembles of non-Hermitian matrices were given by Ginibre |L|. These are the 
ensembles of matrices with arbitrary real, complex, or quaternionic entries. Ginibre gave 
joint probability distributions for the eigenvalues for the complex and quaternionic cases, and 
succeeded in obtaining correlation functions in the complex case, while correlation functions 
for the quaternionic case were found later 0. The purpose of this paper is to extend the 
correlation functions for the quaternionic problem to the weakly non-Hermitian case [^], as 
well as to introduce a fermionic mapping to simplify the computation of these correlation 
functions. Further, the mapping also permits us to derive the 4-point, and higher, correlation 
functions, that were only conjectured before. 

Although Ginibre's ensembles are interesting in themselves, they are also closely con- 
nected with the chiral random matrix ensembles that appear in QCD || and some condensed 
matter systems 0. Knowing the eigenvalue correlation functions in the non-Hermitian en- 
semble, one can easily determine correlation functions in the corresponding chiral ensemble. 
Further, the weakly non-Hermitian versions of these ensembles are of interest in open 
quantum systems; there exists a study using supersymmetric techniques of the eigenvalue 
distribution in the weakly non-Hermitian version of the ensemble considered in this paper 

One interesting feature to observe in the two-level correlation function is the crossover 
from a non-monotonic correlation function with algebraic tails in the limit of very weak 
non-Hermiticity to a monotonically decaying correlation function with Gaussian tails in the 
limit of strong non-Hermiticity. 



1 



Another interesting property that the symplectic non-Hermitian ensemble exhibits is a 
depletion of the eigenvalue density near the real axis; this could be guessed at by looking 
at the joint probability distribution (j.p.d.) derived originally by Ginibre (given in equation 
(§) below). The depletion was also found numerically ||. 

Consider an arbitrary N-by-N matrix of quaternions. This is equivalent to a 2N-by-2N 
matrix M with complex entries. Let M be chosen from an ensemble of such matrices with 
Gaussian weight 

P(M) = e-^ M " M UM (1) 

This defines the strongly non-Hermitian ensemble. The eigenvalues of M come in complex 
conjugate pairs; for every eigenvalue z = x + iy, there is an eigenvalue z = x — iy. Let the 
matrix M have eigenvalues Zi,z iy % = 1...N. Then the j.p.d. of the eigenvalues is given, up 
to a constant factor, by 

' IJe"**!* - Zi\ 2 Y[{zi - Zj)(zi - Zj)(zi - z 3 )(zi - z^JldZidZi (2) 



i<3 



where dzdz = 2dxdy. 

In section II, a fermionic mapping is introduced to write equation (|2]) as a correlation 
function in a fermionic field theory. The mapping is then used to calculate the eigenvalue 
density. In section III, we introduce the weakly non-Hermitian ensemble and calculate 
eigenvalue density for that ensemble. In section IV, multi-eigenvalue correlation functions 
are calculated for both strongly and weakly non-Hermitian cases. The calculations in section 
III and IV are simple extensions of the calculation given in section II. For this reason, the 
calculation in section II is given in the most detail, while the other calculations are sketched. 

For the strongly non-Hermitian case, the main results are equations (^,^) for the eigen- 
value density, and equation ( fHQ for the two-level correlation function. For the weakly non- 
Hermitian case, the main results are equation fl42|) for the eigenvalue density and equations 
(B8|,E9|) for the two-level correlation function. 



II. FERMIONIC MAPPING AND GREEN'S FUNCTION 

In this section we develop the fermionic mapping for the j.p.d. of the strongly non- 
Hermitian ensemble. First, we write equation (0) as a correlation function in a fermionic 
field theory. Then, for convenience we shift to radial coordinates, making a conformal 
transformation. Finally, we integrate over all but one of the Zj to obtain the eigenvalue 
density. The integral over the Zi is done inside the correlation function; only after doing the 
integral is the correlation function evaluated. This amounts to commuting the order of doing 
the integral and evaluating the correlation function, and is the essential trick used in this 
section. In section IV we will demonstrate how to obtain multi-level correlation functions 
by a simple extension of the procedure of this section. 
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A. Fermionic Mapping 

First, let us show that equation (0) can be written as a correlation function in a two- 
dimensional fermionic field theory. A similar fermionic mapping was demonstrated previ- 
ously for the Hermitian orthogonal and symplectic ensembles ||. Let the field ip(x) have 
have the action 

S=-Jdzdz^(z)di)(z) (3) 

Note that we are using only one chirality of fermionic field. Consider a correlation function 
of this field, such as 

27V 

(li^a^ih)) (4) 



This correlation function is equal to 

,2N x ,2N 

(n(oi-ot))(n(6i-6i 



2N 

n (cij - 6, 



(5) 



Let us consider a specific choice of bj, with bj = Le 2m 2N . In the limit L — > oo, we find that 
equation (|5|) reduces to 

2N 

L-™*-" Ilia, -at) (6) 
Comparing this to equation (||D we realize that equation (0) can be written as 

1 2N N 

- hm L 2N2+N ((l[ 4%)) (J! U(zj, ZjMzjWMdZj dzj)) (7) 
00 3=1 j"=i 

where 6j = Le 2m ™ . and £/(z, z) = e~ zz (z — z). When integrating over Zj, the limit L — > oo 
must be taken before doing the integral over Zj. 

Now we will make a conformal transformation to radial coordinates. Write z — e w 
and z = e w . Let w = t + i6. The action for the fermionic field is unchanged under this 
transformation, but we must change equation (^) as the field ip has non-vanishing scaling 
dimension and conformal spin. Equation ([7|) gets replaced by 

i 2N N 

— hm L 2N2 ((Y[ ^( Vj )) (J] &U{e»*, e^w^iw^dwj d Wj )} (8) 

' 3=1 3=1 

where u,- = logfejL + 2m -At. 

Now, we will introduce Fourier transforms for the creation and annihilation operators. 
We will write ip(wi) = J2k ekw * a (k) and ip^ = ^2 k e~ kWi a^(k). In the limit L — > oo, the only 
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states involved in equation @ are those with k = 1/2, 3/2, 5/2, 2N — 1/2. If there are 
excitations in states with higher k, they will vanish in the large L limit. 
Then, we can rewrite equation fl8|), up to factors of order unity, as 



i 27V-1/2 AT _ 

Vf<( II ot W) n(E4K^*^e^(^,^)dWid Wj )) (9) 

JV ' fc=l/2 j=l k,k' 

In equation ([J), consider integrating over Wj,Wj for some given set of j = 1...M. We will 
do the integral inside the correlation function. The integral 

. M _ 

/ II (E a{k)a{k') e w ' k e w ^' e^U{e w ^, e w ')dWj d Wj ) (10) 

J j=l k,k' 

is equal to 

M (11) 

where the operator O is defined by 

= ]Ta(A;)a(£; + l)47r(£; + l/2)! (12) 

k 

Therefore, if we integrate over all eigenvalues in equation (|9|), we obtain 

2JV-1/2 

Z =WA II a\k))0 N ) = U(4n(k + l/2)\) (13) 

JV ' fc=l/2 k 

where the product extends over k = 1/2, 5/2, 9/2, 2N — 3/2. 

B. Eigenvalue Density 

To calculate the density of eigenvalues, we must integrate over all, except one, of the 
coordinate pairs Wi,Wi in equation (H). Using equation ([TT]) , and normalizing with equation 



(Pf), we wish to compute 

__ _ 2JV-1/2 

U(e w ,e w )e t J2 e mw e m ' w (( U at(h)) N " 1 a(m)a(m'))dw dw 

m,m' V k=l/2 ' . . 

(iV — 1)!Z ^ ^ 

It may be verified that the correlation function appearing in the sum of equation (|I~4| ) is 
nonvanishing only if either m — 1/2 is even, m' — 1/2 is odd, and m < m', or if m! — 1/2 
is even, m — 1/2 is odd, and m' < m. In the first case, with m < m r , the contribution to 
equation ([14]) is 

hi R (k+\/2)\ + 1 /2)!e^v / ^(. - z)z m z m 'dwdw (15) 
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where the product over k extends over 

k = 1/2, 5/2, 9/2, ... (16) 
and the product over I extends over 

I = 1/2, 5/2, 9/2. ..m - 2, m + 1, m + 3, m' - 2, m' + 1, m' + 3, 2N - 1/2 (17) 
This is equal to 

1 1 



4tt (m- l/2)!!(m'- 1/2)!! 
In the second case, with m' < m, the result is 
1 1 



e z Vzz(^-z)z m z m d?ydt« 



(18) 



4?r (m- l/2)!!(m'- 1/2)!! 



— z)z m z m dw dw 



(19) 



We can obtain the eigenvalue density p(z, z) by adding equations ( |I8"| , |I9~D and summing 
over m, m'. Shifting m and m! by one-half, and changing from dwdw to dzzw, we find that 
the final result for the eigenvalue density p(z, z)dzdz is 



p(z, z)dzdz = — e zz (z — z)G(z, z)dz dz 



where the Green's function G(z, z) is given by 



G(z, z) 



E 



m<m';m=0,2,4,...;m'=l,3,5,. 



mv.m 



[z m z m -z m z m 



(20) 



(21) 



C. Discussion 



Let us now look at the properties of equation (pH|). We will discuss in turn the normal- 
ization of the density; the way the density depends on x and y separately, where z = x + iy\ 
an integral representation for the density; the circular law; and the depletion of density near 
the real axis. 

First, consider the normalization of the single particle density. It is automatic from the 
above derivation that the eigenvalue density is properly normalized, although one must be 
careful about defining the normalization depending on whether one is counting total number 
of eigenvalues or total number of pairs of eigenvalues. The normalization is defined such 
that / p(z, z)dzdz = N. 

Next, writing z = x + iy and z = x —iy, one can show by differentiating the power series 
in equation (|2~1~1) that d x G(x + iy,x — iy) = ^^-G{x + iy,x — iy) = 2xG(x + iy,x — iy), for 
large N. This implies that 

G(x + iy,x-iy) = e x2 f(y) (22) 
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and therefore p = -^2yie~ y2 f (y) , for some function /, so the interesting properties of the 
eigenvalue density are contained in f(y). Later we will discuss the properties of f(y) for 
small y, and show that there is a depletion of the density of eigenvalues near the real axis. 



Using equations (|2Qj21|j22|) , we can derive an integral representation for p. We can use 
equation (|22| ) to write 

G(z, z) = e zz G(z - z, 0) (23) 
Then equation (plD implies that 

G(?-z,0)= £ \fi-z) m (24) 

m=l,3,5,... m!! 

This is equal to 

oo _ 

r , Sin (i^ t ) e -* 2 /2 dt (25 ) 

J % 



Using this integral representation in equation ( p0[) we find that 

1 °° 

p(z, z) = —2y J Sm(2yt)e- t2 / 2 dt (26) 
o 

Let us now consider the circular law. For large y, equation ( p6|) reduces to 

p(z, z)dzdz — > — dzdz (27) 
Atx 

So, the density tends to a constant for large y. However, this integral representation is 
valid only for iV infinite; for finite N, the density tends to a constant only within a disc of 
radius \/2N, and vanishes outside the disc. This is the well-known circular law The 
vanishing of the density outside the disc is easy to see from the power series representation. 
For finite N the highest power of (zz) appearing in equation (^l|) is roughly 2N and so p 
will be exponentially small for (zz) > 2N. 

Note that the total density in the disc is correct. The area of a disc of radius R is 2ttR 2 , 
where we are using the measure dzdz = 2dxdy. The density is 7-. So, the number of 
particles in a disc of radius V2N is indeed N, as desired. 

For small y, we find that p is reduced below the expected result. Such a reduction was 
found numerically before Q. In the figure, we graph the eigenvalue density as a function of 
y, for x — 0, for a system of 100 particles. 



III. WEAKLY NON-HERMITIAN CASE 

Now we will consider the weakly non-Hermitian version of the ensemble given above. In 
the weakly non-Hermitian random matrix ensemble, we again consider an arbitrary N-by-N 
matrix of quaternions, H, but use a different Gaussian weight. Let H = Hh + H a , where 
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the Hh is Hermitian and H a is anti-Hermitian. Then, we chose the matrix H with Gaussian 
weight 



e -N Tl{H ^ Hh) _N_a Tl{H t aHa) ^ 

where a is some constant. In the large a limit, this reduces to the Gaussian Symplectic 
Ensemble. For finite a, the weight in equation (p8|) is chosen to make sure that the imaginary 
part of the eigenvalues scales with N in the same way as the level spacing. 



If the matrix H is chosen with weight given by equation (p8|) , then the j.p.d. of equation 
gets replaced by 



I JJ e -^?-^|^ JJ^. _ Zj )( Zi -zj){zi - z^Zi-z^Hdzidzi (29) 

i i<j i 



The only difference in the weakly non- Hermitian case is that the Gaussian function of eigen- 
value position e~ ZiZi is replaced by e -3 ^ ~ Na Vi . 

We have not found equation ( P9"D previously in the literature. This equation can be 
derived most easily as follows: write an N-by-N matrix of quaternions, H, as 

H = X~ l TX (30) 

where X is a quaternion matrix such that X^ 1 = X\ and T is an upper triangular matrix 
of quaternions. This procedure is a Schur decomposition, and is possible since the field of 
quaternions, like the field of complex numbers, is algebraically closed. 

The eigenvalues, Zi, can be obtained from the diagonal elements of T; each diago- 
nal element of T is a quaternion, which is associated with a pair of complex conjugate 
eigenvalues z iy Zi. If a given diagonal element of T is Tj = A + Bi + Cj + Dk, then 
Zi = A± WB 2 + C 2 + D 2 . 

The Jacobian associated with this change of variables is n ( z i ~ Zj){ z i ~ Zj). Further, 

i<j 

e -N TT{H lH h )-^Tr(HlH a ) = £ -f TV(T+T h )-^Tt(4T a ) ( 31 ) 

where Th,T a are Hermitian and anti-Hermitian parts of T. The integral over the elements 
of T above the diagonal can be done trivially as this integral is Gaussian. The integral over 
the diagonal elements of T includes a Gaussian factor and a factor from the Jacobian. This 
integral is exactly the integral over the j.p.d. of equation (p9"|) . 

Given equation (|29D , we could follow the procedure of the previous section. However, we 
would run into some difficulties which are purely technical. The problem is that, while in the 
strongly non- Hermitian case the eigenvalue density is independent of x, it is not independent 
of x in the weakly non- Hermitian case. This makes the power series expansion very awkward. 
We will find it convenient to change to a different geometry, given in equation (|32|) below, for 
the weakly non-Hermitian case. Let me again stress that the reason for choosing a different 
geometry is purely technical, to simplify the math. 

An analogous simplification is often used in the Hermitian ensembles. For example, 
consider the Gaussian Symplectic Ensemble in the large iV limit. The eigenvalue density is 
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a function of energy, but if one appropriately scales all energies by the local level spacing, 



it is simpler to obtain correlation functions from the Circular Symplectic Ensemble [JIO 



Let us introduce new coordinates, z = + ir and z = <f> — ir, where is periodic with 
period 2n. Let us replace equation (p9|) by 

1 _ a jvV ( e * z ' — e%Zi ) 2 tt ( e%Zl ~ e tZj )(e tZi — e lZj )(e lZi — e tZj )(e tZi — e tZj ) tt , 

/U| 11 e P -2i<f>i 11 p —2itf>i—2i<f)j \_\_ z i z i \ ) 

■ i i<j i 

This describes a system of N pairs of levels, with average level spacing 2ir/N . The imaginary 
part of the level is of order 1/N, so it is of order the level spacing. For large a this reduces 
to the Circular Symplectic Ensemble. For finite a and large N, we expect that the ensemble 
of equation (|32|) reproduces the behavior of the ensemble of equation ( p9|) within a small 
neighborhood of some given energy, just as the Circular Symplectic Ensemble reproduces 
the results of the Gaussian Symplectic Ensemble within a neighborhood of a given energy. 

The next step is to write equation (^) as a correlation function in a fermionic field 
theory. We will introduce N creation operators at r = +oo and iV creation operators at 
r = — oo. We find that the desired correlation function is 

-i N N 

— Km e^((n^(&i)) (Il^fa)) (n^)^(^)^(%)d^d^)) (33) 
°° i=i i=i i 

where bj = + zL and Cj = ^ — z'L and Z7(r 3 -) = e~ aAr2r j (e rj — e~ rj ). For large N, we can 

write [/(r,-) = e~ aN2r hr r 

Now, we will introduce Fourier modes for the creation and annihilation operators, writing 
V'(^) = J2k elkZa (k) an d ^(z) = J2k e ~ lkZ(l (.k). In the limit L — > oo, the only states involved 
in equation (|33|) are those with k = —N + l/2, —N + 3/2, N — 1/2. If there are excitations 
in states with higher k, they will vanish in the large L limit. Then, equation (|33|) can be 
written as 



N-l/2 

-lim(( I] «t(A;)) n(E«(^)«(^)e^e^f/(r,)d^d^)) (34) 

fc=-A r +l/2 j fc,fc' 



N\ L 



As in the previous section, we will integrate over some set of Zj, for j = 1...M, inside the 
correlation function. The integral 



f ]J(Y, a ( k H k> ) eikZjeik ' ZjU ( r M^j dz j) ( 35 ) 

J j=l k,k> 

is equal to 

OS (36) 

where the operator O w is defined by 

O w = 8(^f /2 Y,( k ^< k H- k )) (37) 
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So, if we integrate over all coordinates z, z in equation (0), we obtain 

1 N-l/2 N-l/2 2 

Z = jn({ n a\k)) (O w f) = U (n^) V2 ke^) (38) 

iv - fc=-jV+l/2 fc=l/2 " 

To obtain the eigenvalue density, we must integrate over all but one of the coordinates in 



equation (|34[) . Using equation (|36|) , and normalizing with equation (j38|) , we obtain 

(39) 



iV-l/2 x 

C/(r) £ e imz e m ' z (( U at(k)) (O w ) N - 1 a(m)a(m'))dz dz 

m,m' V fc=-Ar+l/2 ' 



(N- l)\Z 

The correlation function in equation (p5| ) is non- vanishing only if m = — m'. Equation ( |39| ) 
is equal to 

N-l/2 

U(r) e im ^' z) G w (m)dzdz (40) 

m=-N+l/2 

where 

^ i<^ 3/v * < 41 > 

In the large N limit, we can simplify equation ( flop by introducing scaled coordinates. 
Let us introduce k = m/N and let us also scale z by a factor of N so that now runs from 
to 2nN . Then we can replace the sum by an integral and obtain 

p(0,r)d0dr = ^)- 3/2 re- ar2 G w (z,z)d0dr (42) 

where 

i 

G w (z,z) = I e ik{z - z) e- k2 ' a \dk (43) 

-1 

As in the previous section, the proper normalization of the above result is automatic from 



the derivation. It is possible to show that equation (|42|) is equivalent to equation (26) in the 
limit of very small a. The qualitative feature of a depletion of eigenvalues near the real axis 
is the same for weak and strong non-Hermiticity. Equation ( |4"2]) may be compared to the 
results of the SUSY calculation 0, and found to agree, with some differences in notation 
between the two calculations. 



IV. MULTI-POINT CORRELATION FUNCTIONS 



The calculation of multi-level correlation functions is quite easy. In equations (pj^4t) 



we must integrate over all except for two, three, or more, of the coordinate pairs Wi,Wi. 
Since the system is a non-interacting fermion system, the multi-point Green's functions 
can be expressed very simply in terms of the Green's function (|21|), using Wick's theorem. 
This permits the two-point correlation function to be easily generalized to a multi-point 
correlation function, as conjectured previously [Q. We will not show this in detail, but 
simply sketch the results, first for the strongly non-Hermitian case and then for the weakly 
non- Hermit ian case. 



9 



A. Strongly Non-Hermitian Case 



First let us examine the strongly non- Hermit ian case, generalizing the results of section II. 
Consider the two-level correlation function, the probability to find a pair of levels at position 
z, z given that there is another pair at position z', z' . Then, the two-level correlation function 
is 

(— )V MV (z - z)(z' - z')(G(z, z)G(z', z') - G{z, z')G(z', z) + G(z,z')G(z', z))dzdzdz' dz' 

4:71 ^ ' 

(44) 

This is just an application of Wick's theorem. 

Let us consider the behavior of equation ( [4*4]) in the limit when both z and z are far 
from the real axis. Without loss of generality, assume that Ke(z') > Re(z). Then, use the 
integral representation of the Green's function to rewrite e~ zz ~ z 2 G(z, z')G(zf, z) as 



e 



i,-,'! 2 ( J e (W) V iV2 dt _ I e (W) 2 /2) ^'- Z )V2 _ j e (^) V tV2 dt ) (45) 



We can find a similar representation for G(z,z')G(z',z). Now, in the limit with z and z' 
both far from the real axis, then either Im(z — z') is large or Im(,2 — z') is large. In the first 
case, equation fl4"5|) is exponentially small because of the factor of e~' 2_z ' . In the second 
case, the integrals over t can be performed in this limit, while e^ z ~ z •* I 2 and e^ z _z) I 2 are 

small. The integral, / e^ z ~ z ^e - ' / 2 dt, is equal to -73= , for large Im(V — z); here we rely on 

" 
the fact that Ke(z' — z)>0. 

So, up to exponentially small terms, equation fl4~5D is equal to 



e H ^'' 2 ^— 7=7— (46) 



z — z' z' — z 



11 11 



z—z' z —z z—z z —z' 



Also, in this limit, if equation ( fi6| ) is not exponentially small, then 
Inserting this result, and similar results for G(z,z')G(z', z), back into equation we find 
that, for both z and z' far from the real axis, the two- level correlation function is equal to 

(^) 2 (l - e-\ z ~ z '\ 2 - e-\ z - z '\ 2 )dzdzdz'dz' (47) 

This is essentially the same as the two-level correlation function found in the complex non- 
Hermitian case [0. 

For z and z' near the real axis, I have examined the behavior of equation ([H]) numerically. 
If Im(z) = Im(z'), then the correlation function is a monotonically decaying function of 
Ke(z — z'), with no signs of any oscillation. The correlation function is exponentially small 
if both z and z' are near the real axis. 
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B. Weakly Non-Hermitian Case 



Now let us consider the two-level correlation function in the weakly non- Hermit ian case, 
the probability to find one pair of levels z, z given that there is another pair at z', z'. As 
before, we must integrate over all except for two of the eigenvalue coordinates. Using the 
scaled coordinates, we find that the two-level correlation function is given by 

n \ / 



16 a 



(48) 



As in the previous subsection, this is just an application of Wick's theorem. 

To examine the behavior of the two-level correlation function, let us integrate over r, r' 
in equation Q48I) , to be left with a function of — cf>'. The result is 



i 



l_J_f (k + k'f (fe _ fc0 2 /(2a) 

4tt 2 32tt 2 J kk' K ' 

-i 

In the limit a — > oo, the integral over k, k! in the above expression can be performed to yield 

-2 

Equation fl5"0|) is the known result for the correlation function in the Circular Symplectic 
Ensemble. It is a non-monotonic function, algebraically decaying for large (p. For sufficiently 
small a, equation ( f49| ) will describe a monotonically decaying function of 0' — 0, but for fixed, 
non-vanishing a, the function will always decay algebraically for large 0. 



V. CONCLUSION 

In conclusion, we have given a simple fermionic mapping for determining the correlation 
functions of the non-Hermitian symplectic ensemble. Although the eigenvalue density was 
found previously using SUSY, the present derivation is simpler and can be more easily 
extended to the two-level correlation function. The two-level correlation function in the 
strongly non-Hermitian case was found to be similar to that for the ensemble of arbitrary 
complex matrices. In the weakly non-Hermitian case, the two-level correlation function 
exhibits an interesting crossover as a function of a. 
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